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Abstract:  Using  hypcrspcctral  measurements  made  in  the  field,  we  show 
that  the  cfTeetive  sea-surfaee  refleetanee  p  (defined  as  the  ratio  of  the 
surfacc-rcflcctcd  radiance  at  the  specular  direction  corresponding  to  the 
downwelling  sky  radiance  from  one  direction)  varies  not  only  for  different 
measurement  scans,  but  also  can  differ  by  a  factor  of  8  between  400  nm  and 
800  nm  for  the  same  scan.  This  means  that  the  derived  water-leaving 
radiance  (or  remote-sensing  reflectance)  can  be  highly  inaccurate  if  a 
spectrally  constant  p  value  is  applied  (although  errors  can  be  reduced  by 
carefully  filtering  measured  raw  data).  To  remove  surface-reflected  light  in 
field  measurements  of  remote  sensing  reflectance,  a  spectral  optimization 
approach  was  applied,  with  results  compared  with  those  from  remote- 
sensing  models  and  from  direct  measurements.  The  agreement  from 
different  determinations  suggests  that  reasonable  results  for  remote  sensing 
reflectance  of  clear  blue  water  to  turbid  brown  water  are  obtainable  from 
above-surface  measurements,  even  under  conditions  of  high  waves. 
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1.  Introduction 

The  remote-sensing  reflectance  of  water  (/?„,  sr  *)  is  defined  as  the  ratio  of  the  water-leaving 
spectral  radiance  (Lh,  W  nm'*  sr'*)  to  downwclling  spectral  irradianee  just  above  the 
surface  (£^0^),  W  m'^  nm'*).  R„  (or  L^)  is  the  basis  for  development  of  remote-sensing 
algorithms  as  well  as  for  satellite  sensor  vicarious  calibration  [I  3],  Because  of  various 
technique  limitations  and  the  random  motion  of  the  water  surface,  accurate  determination  of 
Rrs  remains  a  challenge  [2-5],  The  measurement  of  /?„  in  marine  environments  usually 
involves  one  of  these  approaches:  1 )  measure  the  vertical  distributions  of  upwelling  radiance 
(Z,„(z))  and  downwclling  irradianee  {Ej(z))  within  the  water,  and  then  propagate  these 
measurements  upward  across  the  sea  surface  to  calculate  R^  [6];  2)  use  one  sensor  to  measure 
Lu  a  few  centimeters  below  the  surface  and  use  another  sensor  to  measure  £j(0^)  above  the 
surface,  and  then  propagate  £„  across  the  surface  to  calculate  [7];  3)  measure  all  relevant 
quantities  from  an  above-surfaee  platform  [1-5,8-12],  and  then  calculate  Lw  (or  R^,)  by 
removing  surfaee-re fleeted  light  (Lsr).  This  third  approach  is  widely  used  in  the  fteld  and  for 
continuous  measurements  [3,4,11,12],  although  each  approach  has  its  own  advantages  and 
disadvantages  [2,10]. 
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When  measurements  are  made  from  above  the  sea  surfaee  [see  Fig.  1(a)],  the  measured 
signal  is  the  total  upwelling  radianee  (Lr),  which  is  the  sum  of  the  water-leaving  radianee  (Lr) 
and  the  surfaee-refleeted  radianee  (Lsr).  It  is  necessary  to  avoid  viewing  surface  foam,  the 
shadow  of  the  platform  structure,  and  obvious  solar  glint  spots.  Some  surfacc-rcflcetcd  light 
(mostly  from  downwclling  sky  radiance,  but  possibly  including  some  sun  glint)  is  inevitable, 
however.  A  correction  is  therefore  required  to  remove  the  surfaee-refleeted  light  from  Lt  in 
order  to  compute  Lw  and  Rrs  [2,4,8].  One  approach  to  the  removal  of  surfaee-refleeted 
radianee  was  proposed  by  Mobley  [13]  (a  similar  description  ean  be  found  in  Morel  [1]).  In 
this  technique,  all  Lsr  is  expressed  as  the  product  of  p  -  an  effective  surfaee  rcfleetanee  and 
sky  radianee  {Lsky)  measured  for  an  angle  reciprocal  to  the  measurement  of  Lr  (see  Fig.  1  in 
Ref.  [13]).  The  value  of  p  depends  on  sea  state,  sky  conditions,  and  viewing  geometry  [13,14]. 
Two  approaches  have  then  been  proposed  for  the  determination  of  p:  One  is  to  derive  the 
value  of  p  from  measured  Lt  and  Lsky  by  assuming  L^  approaches  0  at  near-infrared 
wavelengths  (e.g.  at  780  nm)  [  1  ];  the  other  is  to  use  a  table  of  p  values  derived  from  numerical 
simulations  with  various  wind  speeds  and  viewing  geometries  [13].  Both  approaches  [1,13], 
however,  assume  that  the  p  value  is  spectrally  constant.  To  minimize  the  impact  of  sun  glint 
on  the  derivation  of  Lw  (or  /?„),  Flooker  et  al.  [2]  and  Zibordi  et  al.  [4]  suggested  filtering  out 
the  higher  measured  total  radiance  {Lt)  values,  and  reasonably  good  results  were  achieved  for 
Lw  in  the  412-555  nm  range  (larger  uncertainties  were  found  at  670  nm  [4]).  Here,  after 
describing  the  general  dependence  of  p,  we  show  with  hypcrspcetral  measurements  that  p  in 
general  varies  with  wavelength,  and  that  the  spectral  variation  ean  be  significant.  We  further 
compare  two  physical-mathematical  approaches  and  a  direct  measurement  scheme  for  the 
removal  of  Lsr  in  deriving  /?„. 


Fig.  1.  (a)  Schematic  illustration  of  abovc-surfacc  measurement  of  Lt.  (b)  Example  of 
roughened  sea  surface  when  looking  down  from  an  above-surface  platform.  The  different 
shades  of  blue  result  from  light  reflected  from  different  parts  of  the  sky. 


2.  Theoretical  background 

When  a  radiance  instrument  takes  measurements  of  spectral  upwelling  radiance  (^/{A.))  from 
an  above-surfaec  platform,  it  collects  not  only  the  radiance  emerging  from  below  the  water 
surfaee  (the  so-ealled  water  leaving  radianee,  Z.h(^)),  but  also  surfaee-refleeted  light  (LsRi)^)). 
If  the  measurement  angle  is  0  from  nadir  and  (p  (azimuth)  from  the  solar  plane  [sec  Fig.  1(a)], 
then  for  a  level  sea  surfaee,  Ls/^)^)  comes  from  the  zenith  angle  0’  =  0  and  the  same  azimuthal 
angle  ((p).  For  the  more  common  situation  of  a  constantly  moving,  roughened,  surface  (sec 
Fig.  1(b)  for  an  example),  and  for  typical  instrument  integration  times  of  order  of  one  second 
or  longer  (integration  time  is  much  shorter  for  multiband  sensors  [3]),  LsRi'k)  actually  eomes 
from  a  large  portion  of  the  sky  (see  Fig.  1  and  Fig.  2  of  Mobley  [13])  and  may  include  solar 
radianee  (sun  glint).  So,  in  general,  the  spectral  upwelling  radiance  measured  from  an  angular 
geometry  (0,(p)  is 
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L,{X,e.(p)  =  \<p. \ 


(1) 


Here  subseript  “f’  represents  the  small  wave  faeet  viewed  by  the  sensor;  Wi  is  the  relative 

weighting  of  solid  angle  of  the  faeet  to  the  sensor’s  field-of-view  solid  angle;  F  is  the 

Fresnel  refleetanee  of  the  faeet;  and  Lsky  is  the  down  welling  radiance  ineident  onto  the 
faeet,  whieh  is  refleeted  into  senor’s  viewing  angle. 

Because  LsR(k)  is  assembled  in  an  unknown  manner  according  to  Eq.  (1),  removal  of 
LsrO^)  beeomes  a  challenge  in  the  field  when  measurements  are  taken  from  above  the  sea 
surfaee  (or  sea-surfaee  remote  sensing  in  analogy  to  satellite  remote  sensing).  For  this 
removal,  to  a  first  order  approximation,  Eq.  (1)  is  simplified  as  [1,13] 

Lr(Ji,0,<p)  =  L„.(Ji,0,<p)  +  p(0,<p)Lsi„XJi,0',<p).  (2) 


Here  LskyiO\(p)  is  the  sky  radianee  in  the  same  plane  as  that  of  Lr,  but  with  0’  the  reeiproeal 
(speeular)  angle  of  0  [2,4, 8,9].  p(0,(p)  is  the  effective  surface  reflectance  that  accounts  for 
reflected  sky  light  from  all  directions  for  the  given  sensor  direetion,  and  is  assumed  to  be 
independent  of  wavelength.  p(0,(p)  equals  the  Fresnel  refleetanee  of  the  sea  surface  only  if  the 
surfaee  is  flat  (without  waves).  Values  of  p(0,(p)  for  various  viewing  directions,  sun  zenith 
angles,  and  wind  speeds  were  evaluated  with  numerieal  simulations  in  [13].  Based  on  these 
simulations,  it  was  suggested  to  use  0  =  40®  from  the  nadir  and  (p  =  135®  from  the  sun  to 
minimize  Lsr  when  measuring  Rrs  in  the  field. 

Comparing  Eqs.  (1)  and  (2)  gives 


p{0,<p)  = 


S  '■  ’’  V' )  Av*. ’’  V, ') 


(3a) 


or. 


_  Lr{A.,0,<p)-LJA.,0,<p) 
L^iX,0\<p) 


(3b) 


Because  Lsky^  in  general  has  different  speetral  shapes  for  different  directions  [1]  (e.g.,  for  a 
clear  sky  at  noon,  Lsky  from  the  horizon  appears  whiter  than  that  from  the  zenith),  p  in  Eq.  (2) 
or  Eq.  (3a)  will  in  general  be  speetral ly  dependent,  especially  when  solar  light  is  inevitably 
reflected  into  sensor’s  viewing  angle  by  roughened  surfaee,  unless  the  sky  is  completely 
overcast. 

3.  Data  and  methods 

To  demonstrate  the  spectral  variation  of  p,  hyperspectral  measurements  over  clear  oceanic 
waters  were  utilized,  where  the  contribution  of  water  to  Lt  is  negligible  in  the  longer 
wavelengths.  The  measurements  were  made  on  Feb.  23,  1997  around  12:50  pm  (local  time), 
for  waters  near  Hawaii  at  21.33  N,  158.16  W.  The  sky  was  clear  with  no  clouds,  the  water 
appeared  blue,  the  wind  speed  was  around  8ms',  and  the  surface  wave  amplitude  was  -2  to 
3  feet  (-“1  m). 

Upwelling  total  radiance  (L/-,  9  scans),  downwelling  sky  radiance  {Lsk^^  5  scans),  and 
“gray-card"  radiance  {Lc,  3  scans)  reflected  from  a  standard  diffuse  reflector  (Spectra Ion®) 
were  measured  with  a  handheld  speetroradiometer  (SPECTRIX  [15]),  which  covers  a  speetral 
range  -360  -  900  nm  with  a  speetral  resolution  -2  nm  and  has  an  integration  time  about  1 .5 
seconds  for  the  collection  of  Lj.  The  orientation  to  measure  Lj  was  30®  from  nadir  and  90® 
from  the  solar  plane.  Lsky  was  measured  in  the  same  plane  as  Lr,  but  at  an  angle  30®  from 
zenith.  Downwelling  irradianec  was  determined  by  assuming  that  the  Spectralon  is  a 
lambertian  reflector,  so  that  Ej  =  nLc/Rcy  with  Lc  the  average  of  the  three  scans  and  Ra  the 
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reflectance  of  the  diffuse  reflector  (--10%).  The  measurement  was  taken  at  the  bow  of  a  large 
ship  with  a  sensor  to  water-surface  distance  about  5  meters.  The  SPECTRIX  has  a  10®  field  of 
view,  which  then  observed  a  surface  area  of--!  m^  for  this  setup. 

To  evaluate  the  value  and  variations  of  the  effective  surface  reflectance  (p),  Eq.  (2)  is 
converted  to  reflectances,  where  the  total  remote-sensing  reflectance  ratio  of  Lt  to  EJ) 
and  sky  remote-sensing  reflectance  (5^,,  ratio  of  Lsky  to  E^)  were  calculated  for  each  Lt  and 
Lsky  scan,  respectively.  From  Eq.  (2),  these  /?„,  and  are  related  as 

R,M>0.(p)^  p(0.(p)SM,0\(p\  (4a) 


or 


K  <P)  =  <P)  -  ^  <py 


Further,  p  is  calculated  as 


p(0.<p)^ 


(4b) 

(5) 


For  the  calculation  of  p,  /?„(X,30®,90°)  (assumed  equal  to  ^„(X,0®,0®),  and  for  X  in  a  range  of 
400  800  nm)  was  estimated  with  the  bio-optical  model  of  Morel  and  Maritorena  [16]  and 

using  a  chlorophyll-a  concentration  of  Chi  =  0.05  mg  m“^,  which  is  an  estimate  based  on 
observations  of  MODIS  for  these  waters  in  February.  The  model  coefficients  of  Morel  and 
Maritorena  [16]  cover  wavelengths  up  to  700  nm.  For  the  study  here,  the  model  coefficients  of 
wavelengths  greater  than  700  nm  arc  considered  the  same  as  that  at  700  nm,  except  for  the 
attenuation  coefficient  of  pure  water,  which  was  replaced  with  the  absorption  coefficient  of 
clear  natural  water  [17]. 

4.  Results 

4.  /  Variation  of  p 

For  this  station,  and  Sr,  arc  presented  in  Figs.  2(a)  and  2(b),  respectively.  Because  the  sea 
surface  is  roughened  by  waves,  as  commonly  encountered  in  the  field,  we  did  not  get  identieal 
Trs  for  the  9  independent  measurements  of  Lj.  This  is  because  each  Lj  measurement  observed 
a  different  sea  surface,  hence  a  different  sky  coverage,  and  thus  a  different  Lsr.  Nor  did  we  get 
identical  measurements  of  the  5  because  the  boat  was  also  constantly  moving,  and  thus  the 
sensor  could  not  maintain  the  exactly  same  angular  geometry  for  the  different  sky-viewing 
measurements. 


Fig.  2.  Measured  TV,  (a)  and  (b). 
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Fig.  3.  p  values  calculated  from  measured  Tn  and  5r,.  Rn  was  modeled  with  Chi  0.05  mg  m  ' 
based  on  the  bio-optieal  model  of  Morel  and  Maritorena  [16]. 

For  illustration  purposes.  Fig.  3  shows  values  of  p  calculated  for  the  9  scans  and  with 
Srs  from  the  first  measurement  used  as  the  denominator  in  Eq.  (5).  It  is  seen,  not  surprisingly, 
that  the  p  values  differ  among  the  different  Lt  measurements.  More  importantly,  the  p  values 
differ  spectrally,  and  this  difference  can  be  as  high  as  a  factor  of  eight  between  400  nm  and 
800  nm  (a  factor  of  five  between  400  nm  and  700  nm).  The  increase  of  p  with  wavelength 
occurs  mainly  because  (1)  collects  Lsr  from  all  directions,  including  the  sun  and  near¬ 
horizon  directions  [recall  the  whitish  patches  in  Fig.  1(b)].  Compared  to  sky  light  from  zenith, 
radiances  from  these  directions  arc  richer  in  the  longer  wavelengths.  (2)  is  measured  from 
one  fixed  angular  geometry,  and  this  lS’„  is  usually  blue  rich  (dominated  by  contributions  from 
Rayleigh  scattering). 


Wavelength  [nm] 

Fig.  4.  Similar  as  Fig.  3,  bul  wilh  two  different  Chi  values.  Green:  Chi  -  0.05  mg  m  blue.  Chi 
-  0.1  mg  m  \ 

To  evaluate  the  impacts  of  incorrect  /?„,  which  were  estimated  from  a  spectral  model  with 
roughly  estimated  chlorophyll  concentration,  on  the  calculated  p  values.  Fig.  4  compares  the  p 
values  calculated  from  the  9  measurements  and  the  first  but  with  Chi  =  0.05  and  0.1 

mg  m  respectively.  For  wavelengths  in  the  range  of  -400  500  nm,  because  /?„  makes 
strong  contributions  to  7V,,  wide  variations  of  p  values  were  found,  which  highlights  the 
limitation  of  calculating  the  effective  p  from  field  measurements  when  the  water  contribution 
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is  high.  For  wavelengths  longer  than  '-SSO  nm,  however,  it  is  found  that  the  impact  of 
different  Chi  values  (then  different  /?„)  on  p  is  nearly  negligible.  This  is  beeause  for  sueh  elear 
waters  phytoplankton  eontribution  to  /?„  is  nearly  negligible  at  the  longer  wavelengths.  This  is 
further  illustrated  in  Fig.  5  via  seatter  plots  between  p(Chl  =  0.05)  and  p(Chl  =  0.025),  and 
between  p(Chl  =  0.05)  and  p(Chl  =0.1).  This  figure  shows  that  Chi  has  very  little  impact  on  p 
values  of  p  >  0.07  (corresponding  to  --5 50  nm  for  the  measurements  in  this  study).  The  same 
results  were  found  when  the  first  was  replaced  by  any  of  the  other  incasurcments  of  Lsky 
(results  not  shown  here). 


P 

Fig.  5.  Scatter  plot  between  p(ChI  ~  0.05)  and  p<Chl  =  0.025),  blue  symbol,  and  bclwccn  (XChl 
=  0.05)  and  p(Chl  =  0.1),  green  symbol. 

If  there  arc  clouds  in  the  sky  (assuming  the  sun  itself  is  not  blocked  by  clouds),  this  p 
value  could  vary  widely  with  wavelength,  because  could  be  measured  from  a  small  portion 
of  the  clear  sky  (very  blue)  or  aimed  at  a  cloud,  while  Lsr  will  include  radiance  from  clouds 
(nearly  white)  and  the  background  blue  sky.  These  results  indicate  that  applying  a  p  value 
calculated  in  the  near  infrared  (c.g.  780  nm)  to  the  shorter  wavelengths  will  cause  large 
uncertainties  in  in  the  blue  bands  [2],  unless  the  measurements  arc  made  under  nearly  ideal 
conditions  (no  clouds,  low  wind,  no  foam  on  surface,  and  very  short  integration  time). 

4.2  Removal  of  Lsr 

The  above  analysis  indicates  that  when  the  sea  surface  is  not  flat,  1)  p  is  not  a  constant  among 
measurement  scans;  and  2)  p  values  change  with  wavelength,  at  least  for  the  longer 
wavelengths  (>  -550  nm)  in  this  study.  With  such  an  observation,  even  if  wind  speed  and 
angular  geometry  arc  all  known  exactly  (note  that  the  effective  p  also  depends  on  the 
orientation  of  waves),  it  will  still  be  a  daunting  challenge  to  accurately  remove  Lsr  via  Eq.  (2) 
or  Eq.  (4).  Earlier,  Hooker  et  al.  [2]  and  Zibordi  et  al.  [4]  proposed  to  filter  out  the  higher  Lj 
measurements  before  applying  Eq.  (4b)  for  the  removal  of  Lsr^  This  technique  is  generally 
supported  by  the  results  showTi  in  Fig.  3,  where  higher  spectral  contrast  of  p  is  found  for  the 
high  p(800)  value  (high  Ly).  However,  because  it  can  never  be  known  exactly  which  Lsky  is 
reflected  into  the  view  of  an  Lj  measurement,  it  is  unclear  how  to  select  a  proper  p  value  that 
is  relevant  for  the  smaller  Lt  measurements,  as  using  the  smallest  Lj  to  derive  L»  via  Eq.  (4b) 
may  result  in  underestimation  of  Lh-  [  1 8]. 

For  this  station,  the  wind  speed  was  about  8  m  s  *,  so  a  p  value  of  0.05  was  assumed  for 
the  angular  geometry  (based  on  Fig.  8  in  Mobley  [13])  and  applied  for  the  calculation  of  R^s 
via  Eq.  (4).  Since  there  were  9  measurements  of  and  5  measurements  of  45  /?„  were 
derived.  Figure  6  shows  the  average  /?„  with  ±  1  standard  deviation  as  computed  from  the  45 
spectra.  As  a  qualitative  check,  the  modeled  /?„  for  Chi  -  0.05  mg  m  ^  is  also  included  in 
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Fig.  6.  It  is  found  that  the  average  from  measurements  match  modeled  /?„  reasonably  well 
for  the  ~400-550  nm  range,  but  there  are  significant  differences  for  the  longer  wavelengths. 
Since  there  are  large  uncertainties  in  the  modeled  Rrs  (resulted  from,  likely,  both  inaccurate 
Chi  value  and  imprecise  Rrs  model),  we  are  not  expecting  the  two  matching  each  other 
exactly.  However,  because  the  water-leaving  radiance  (or  /?„)  of  such  waters  is  nearly 
negligible  at  longer  wavelengths,  it  ean  be  safely  argued  that  the  /?„  derived  from  Eq.  (4)  is 
overestimated  for  those  wavelengths.  This  observation  is  consistent  with  Fig.  13  (left)  of 
Mobley  [13]. 


Fig.  6.  Average  (solid  blue)  and  the  one  standard  deviation  (dotted  lines),  calculated  using 
Eq.  (4b).  Green  line  is  modeled  R„  with  Chi  =  0.05  mg  m  ’  using  Ihc  Morel-Marilorcna  hio- 
oplieal  model  [16]. 

The  commonly  measured  properties  (except  grey  card)  are  Lt  at  angle  (0,cp)  and  Lskv  at 
(0\(p).  Also,  because  the  actual  total  Lsr  is  not  measured  directly,  wc  re-write  Eq.  ( I )  as 

L,(X,0,<p)  =  L^(A.(9,(3)+  F(0,(p)  L^.(X,0\<p)  +  Y.  w,F(0, \(p, \0,(p)  (A.O, \<p, '). 

j-l 

(6) 

Here  wo  represents  the  weighting  of  sky  light  coming  from  the  specular  direction  (0\(p),  and 
the  sum  now  includes  sky  light  from  all  other  directions.  Further,  since  sky  light  from  the 
specular  direction  (0’,(p)  dominates  Lsr  from  the  fleld-of-view  centered  at  (0,(p)  [13],  Eq.  (6)  is 
approximated  (by  setting  wo  =  1 )  as, 

Lj.{A,0,<p)x  L„  {A,0,(p)  +  F(0,(p)  L^XA,0\(p)  +  '^  F{0,\(p/,0.(p)  L^^  (A,0,\(p,\ 

i  I 

(7) 

or,  in  terms  of  reflectance, 

T^{A,0,<p)  ^  R^(A,0,<p)-\~  F(0,(p)  u;  F{0^  \<p^  \  0,(p)  \  (p^ "). 

j  > 

(8) 

Now  in  Eq.  (8)  both  T„  and  for  the  specular  direction  are  directly  detenmned  from 
measurements.  F(0,(p)  is  the  Fresnel  reflectance  of  water  surface  for  (0,(p),  which  is  known  for 
a  given  angular  geometry.  For  the  calculation  of  Rrs  from  Eq.  (8)  for  any  measurement  of  Lt 
and  Lsk\^  it  is  thus  necessary  to  determine  the  last  term  on  the  right-hand  side  of  Eq.  (8). 
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Because  it  is  not  known  yet  how  this  last  term  varies  spectrally,  this  tenn  is  assumed  for 
expediency  to  be  spectrally  independent  [9].  Thus  Eq.  (8)  becomes 

TrM^O,(p)  ^  R„(X,0.(p)’¥  F(0,(p)  S,M>0\(p)^^(0.(p\  (9a) 


or  [9], 


RM^O,(p)^T^MM-nO.(p)S^(?i,0\(p)-^(0.(p).  (9b) 


Thus,  for  each  set  of  spectral  Trs  and  there  is  a  spectrally  constant  value  (A,  a  bias)  that 
must  be  determined  before  Rr^  ean  be  derived.  For  oceanic  waters  where  R„  is  negligible  in 
the  red  and  near  infrared,  A  ean  be  estimated  by  assuming  Rr^  near  750  nm  is  0  [19].  For 
coastal  turbid  waters,  however,  this  assumption  is  no  longer  valid.  For  such  environments,  one 
approach  [19,20]  is  to  model  the  spectral  R„  as  a  function  of  spectral  inherent  optical 
properties  (lOPs),  and  then  solve  for  A  by  comparing  modeled  spectral  R„  with  spectral  /?„ 
derived  from  Eq.  (9b)  using  all  measured  spectral  information  (so-ealled  spectral 
optimization)  [21-24]. 

Basically,  for  optically  deep  waters,  the  spectral  /?„  can  be  conceptually  summarized  as 
Rr,{^,0,<p)^  Fun(a(X), b^(X),0,<p\  (10) 


where  a(X)  is  the  absorption  coefficient,  and  hf,(X)  is  the  baekscattering  eoeffieient.  The 
inherent  optical  properties  a(X)  and  hff(X)  can  be  modeled  with  bio-optieal  models  of  optically 
active  components  [22,24,25],  so  that  Eq.  (10)  becomes  explicit  functions  as 

RJA,^ 0. (p)  ^  Fun{a^.(A, ), ),  P, C,  X, P, (p\ 

Funia^  C,  ^  ^  ^  ^ 

RM„M^  I  0. (p). 


Here  X\  to  Xn  are  the  sensor’s  wavelengths,  and  are  the  knowTi  absorption  and 
baekscattering  coefficients  of  pure  seawater,  and  P,  G,  and  X  represent  the  magnitude  of  the 
absorption  eoeffieient  of  phytoplankton,  gelbstoff,  and  the  baekscattering  eoeffieient  of 
particles,  respectively.  To  derive  the  value  of  A  in  Eq.  (9b),  an  objective  function  is  defined  as 


Err  = 


^  I  iKOO  ^ 


(12) 


with  R^  from  Eq.  (11)  while  /?„  from  Eq.  (9b).  represents  the  average  of  an  array 

between  400  nm  and  675  nm.  The  upper  bound  of  wavelength  (800  nm)  can  be  extended  to  a 
longer  wavelength  for  turbid  lake  or  river  waters  when  sensor  has  measurements  in  those 
wavelength  ranges.  Err  is  then  a  function  of  4  variables  (P,  G,  X,  and  A)  for  optically  deep 
waters,  and  they  are  derived  numerically  when  Err  reaches  a  minimum  spectral  optimization 
[22,26].  Rrs  is  therefore  computed  by  applying  this  numerically  derived  A  to  Eq.  (9b).  Note 
that  in  the  correction  of  Lsr  the  foeus  is  the  estimation  of  A,  although  values  of  P.  G.  and  X  are 
also  determined. 

For  the  measurements  at  this  station,  again,  45  spectral  P„  were  detennined  with  this 
spectral  optimization  method,  and  their  average  and  standard  deviation  are  presented  in  Fig.  7. 
The  overestimations  of  R„  in  the  longer  wavelengths  (>  --550  nm)  are  generally  removed,  as 
compared  to  Fig.  6.  At  the  same  time,  the  average  matches  the  modeled  (with  Chi  = 
0.05  mg  m”^)  very  well  in  the  -^00-550  nm  range,  although  it  is  not  our  intension  (the 
measured  and  Chl-modeled  R„  do  not  necessarily  represent  the  same  water  environments). 
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Fig.  7.  Similar  as  Fig.  6,  but  /?„  was  calculated  based  on  Eq.  (9)  and  with  a  spectral 
optimization  scheme.  Green  line  is  modeled  Rn  with  Chi  “  0.05  mg  m  ^  using  the  Morcl- 
Maritorena  [16]  bio-optical  model. 

To  further  lest  the  above  evaluation  and  the  optimization  approach  of  removing  Lsr^  new 
measurements  (September  13,  2010;  ~1 1  am  local  time)  were  carried  out  (with  SPECTRIX) 
over  turbid  river  water  (Pearl  River,  Mississippi,  USA.  Figure  8  shows  color  photos  of  the 
water  and  sky  when  measurements  were  taken).  This  shallow  (-0.5  m)  and  very  turbid  water 
makes  it  nearly  impossible  to  obtain  /?„  from  measurements  of  vertical  profiles  of  Ly  and 
[6].  During  the  experiment,  the  surface  was  calm  [sec  Fig.  8(a)]  and  the  sky  was  blue 
[Fig.  8(b)]  with  some  thin  cirrus  clouds.  Two  different  measurement  schemes  were  carried 
out.  One  followed  the  traditional  scheme  [10]  that  measures  Lj  and  Lskv  (see  Section  3), 
with  0  =  30°  from  nadir  and  (p  =  90°  from  the  solar  plane,  and  the  sensor  to  water-surface 
distance  was  -1  meter  (the  sensor  then  covered  a  surface  area  -0.05  m").  were  derived, 
separately,  from  these  measurements  following  the  simple  approach  [Eq.  (4b),  p  =  0.022  is 
used  for  calm  surface.  Rrs'Simp  in  the  following]  and  following  the  optimization  approach 
(Rrropt  in  the  following)  mentioned  above. 


Fig.  8.  Color  photos  of  the  river  water  (a)  and  sky  (b)  measured  on  September  13.  2010,  -  I  I 
am  local  time. 

Another  measurement  followed  a  novel  scheme  proposed  by  Ahn  et  al  [27],  where  a  small 
black  tube  (-4  cm  in  diameter)  was  placed  in  front  of  the  sensor  to  block  Lsr  (see  Fig.  9  for  a 
schematic  illustration).  When  L^was  measured  the  tube  was  dipped  just  below  the  sea  surface 
(-5  cm)  while  the  sensor  itself  was  kept  above  the  surface.  Therefore  there  will  be  no  Lsr  into 
the  sensor  in  this  setup  and  the  instrument  records  a  direct  measurement  of  L^.  Rrx  (Rrs^direct 
in  the  following)  was  then  derived  as  the  ratio  of  measured  Lh  to  (from  measurement  of 
LgI 
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sensor 


Fig.  9.  Scheme  to  measure  Lw  directly  (re-drawn  from  Ahn  cl  a)  [27]).  The  open  hox  is  a  black 
tube  to  block  surface  reflected  light,  which  is  inserted  just  below  (-5  cm)  the  surface  when 
measuring  Lw. 

Figure  10  shows  the  derived  Rrs  irom  the  three  measurement-determination  sehemes;  blue 
is  Rrs-directy  green  is  Rr^-opt^  and  eyan  is  Rr,‘Simp.  The  three  Rr,  eurves  show  similar  speetral 
shapes,  whieh  are  typieal  of  turbid,  high-CDOM  river  waters  (note  the  yellow-brown  eolor  in 
Fig.  8).  Rrs’Simp  is  eonsiderably  higher  than  both  Rr^-direct  and  Rr^-opt,  suggesting  ineomplete 
removal  of  Lsr  even  for  this  quite  ealm  situation  (it  may  be  that  some  sun  glitter  could  not  be 
eompletely  avoided  for  the  (30°,90°)  viewing  geometry  and  integration  times  of  -1-2 
seeonds).  On  the  other  hand,  Rrs-direct  and  R^-opt  are  very  eonsistent  aeross  the  --400-850  nm 
range,  with  a  eoeffieient  of  variation  about  -11%  (whieh  is  about  46%  between  and 

Rrs-direct).  The  slight  negative  /?„  (both  R^rdirect  and  Rr^-opt)  for  wavelengths  shorter  than 
400  nm  may  result  from  a  eombination  of  1)  SPECTRIX  has  lower  signal-to-noisc  ratio  for 
wavelengths  shorter  than  —400  nm  [15],  and  2)  the  extremely  low  upwelltng  signal  in  the 
blue-to-UV  wavelengths  of  this  CDOM-rieh  water.  Nevertheless,  the  deduced  Rrs  of  this 
turbid  water  (along  with  the  result  of  blue  oeeanie  waters)  strongly  indicates  that  Eq.  (9b)  with 
an  optimization  seheme  to  determine  the  value  of  A  is  adequate  in  obtaining  Rrs  in  the  field 
when  measurements  are  made  above  the  sea  surfaee  under  un-ideal  eonditions  and  that  Lsr  is 
not  bloeked  during  measurements.  However,  neither  Rrs-direct  nor  Rrs-opt  are  error  free, 
beeause  Rrs-direct  eneounters  some  self-shading  and/or  eontributions  from  refleetanee  inside 
the  tube,  while  Rrs-opt  suffers  from  the  approximation  from  Eq.  (6)  to  Eqs.  (7)  and  (9a). 


400  500  600  700  800 

Wavelength  [nm] 


Fig.  10.  Comparison  between  direclly  measured  R„  (blue  line)  of  water  showing  in  Fig.  X(a) 
and  /?„  obtained  after  couecling  surface- reflected  light. 
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5.  Conclusions 


Using  measurements  from  a  elear-water  station,  we  demonstrated  that  the  effeetive  surfaee 
refleetanee  (p)  varies  not  only  with  each  measurement  scan  but  also  with  wavelength. 
Consequently,  applieation  of  a  speetrally  constant  p  value  for  the  removal  of  Lsr  from  above- 
surface  measurements  is  a  crude  approximation,  especially  if  the  sea  surface  is  significantly 
roughened  by  waves  and  the  sensor  has  a  long  integration  time  (as  do  most  high-spectral 
resolution  sensors).  Earlier  studies  [2,4,18]  have  shown  that  it  is  wise  to  filter  out  the  higher 
Lf  measurements  before  the  derivation  of  Lk  when  the  simple  formula  [e.g.,  Eq.  (4b)]  is  used 
for  the  derivation.  Here  we  show  that  for  clear  to  turbid  waters,  a  spectral  optimization 
scheme  [28]  is  also  adequate  to  remove  Lsr  \i\  Lj  measurements  and  derive  reasonable  /?„. 
Further,  the  scheme  to  block  Lsr  by  equipping  a  black  tube  in  front  of  the  sensor  and  dipping 
it  just  below  the  surface  shows  promise  to  obtain  reliable  measurement  of  Lw  without  the 
difficult  post-processing.  Further  effort  is  required  by  the  remote-sensing  community  to 
evaluate  these  approaches  for  a  wide  range  of  environments  and  measurement  conditions  (e.g. 
Hooker  ct  al.  [2]  and  Toole  et  al  [7])  and  then  to  establish  a  consensus  for  the  optimum  way  to 
determine  /?„  in  the  field  when  measurements  are  made  from  an  above-surface  platfonn, 
especially  for  situations  such  as  turbid  waters  and  partly  cloudy  skies. 
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